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^ ■ Abstract 



The electroweak contributions to the forward-backward asymmetry in the production of top- 
l-- quark pairs at the Tevatron are evaluated at 0{a^) and 0{aaf). We perform a detailed 

analysis of all partonic channels that produce an asymmetry and combine them with the QCD 
contributions. They provide a non-negligible fraction of the QCD-induced asymmetry with the 
same overall sign, thus enlarging the Standard Model prediction and diminishing the observed 
deviation. For the observed mass-dependent forward-backward asymmetry a 3a deviation still 
tv^ ' remains at an invariant-mass cut of M^t > 450 GeV. 
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1 Introduction 

The investigation of the top quark at the Tevatron has substantially contributed to precision tests of 

QCD and the electroweak theory. Besides the valuable set of top-quark observables like mass, width, 

cross section, which are fully consistent with the SM, the measured forward-backward asymmetry 

ApB of top-pair production [TJ[2] is larger than expected from the Standard Model (SM) prediction 

and has led to speculations on the presence of possible new physics. 

Two options for the forward-backward asymmetry have been used in the experimental analysis, 

with the definitions 

<7{Ay > 0) - a{Ay < 0) 



and 



"^^^ a{Ay > 0) + <T{Ay < 0) ^^^ 

App _ <^{yt > 0) - (TJyt < 0) ,„. 

^^^~a(y,>0) + a(y, <0) ^'> 

given in [3] reporting the recent CDF result. Ay is defined as the difference between the rapidity yt 
and yt of t and i where the proton direction defines the beam axis. Ay (not yt) is invariant under 
a boost along the beam axis, thus it is the same in the partonic and in the hadronic rest frame. 
The recent values for the inclusive asymmetry obtained by CDF [3| are 

A^^B =0.158 ±0.075, (3) 

A% = 0.150 ±0.055. 

The LO predictions of the Standard Model originate from NLO QCD contributions to the differential 
cross section for tt production that are antisymmetric under charge conjugation [HIS], yielding values 
for Apg{A^g) around 7%(5%) (see e.g. [6]). The observed difference between the measurement 
and the prediction has inspired quite a number of theoretical papers proposing various new physics 
mechanisms as potential additional sources for the forward-backward asymmetry (sec for example [71 
[5] and references therein). 

The importance of identifying signals from possible new physics requires a thorough discussion of 
the SM prediction and the corresponding uncertainty. At present, the theoretical accuracy is limited 
by the missing NNLO contribution from QCD to the antisymmetric part of the ti production cross 
section. Besides the strong interaction, the electroweak interaction gives rise to further contributions 
to the tt forward-backward asymmetry, through photon and Z exchange at the tree level as well 
as through interference between QCD and electroweak amplitudes at NLO in both interactions. 
Although smaller in size, they are not negligible, and a careful investigation is an essential ingredient 
for an improved theoretical prediction. 

In this paper we perform a detailed analysis of the electroweak contributions to the forward- 
backward asymmetry in tt production based on the evaluation of all partonic channels that produce 
an asymmetry both at the tree level and at NLO, and combine them with the QCD contributions. 
We apply the calculation to both types of asymmetries given above in ^ and ([2]). Moreover, we 
present results for A"^ also with a cut Mtt > 450 GcV on the invariant tt mass as well as with a 
rapidity cut | Ay | > 1 , for comparison with the experimental values given in [3] , 

ApsiMti > 450 GcV) = 0.475 ± 0.114, A^lgHAyl > 1) = 0.611 ± 0.256, (4) 

where in particular the result for the high invariant-mass cut exhibits the largest deviation from 
the QCD prediction. 



2 Calculational basis 

At leading order the production of ti pairs in pp collisions originates, via the strong interaction, 
from the partonic processes qq —j' tt and gg -^ ti, which yield the 0{al) of the (integrated) cross 
section, i.e. the denominator of ApB in ([1]) and ^. The antisymmetric cross section, the numerator 
of ApB, starts at Oia'l) and gets contributions from qq — >■ ti{g) with q = u,d (the processes from 
other quark species, after convolution with the parton distributions and summation, are symmetric 
under yt — > —yt and thus do not contribute to Apb) as well as from qg — )■ ttq and qg -^ ttq. 

Writing the numerator and the denominator of ApB (for either of the definitions ([T|) and ([2])) 
in powers of Us we obtain 

(5) 






D ajDo + alDi 



= 7f (VVi + as{N2 - NiDi/Do)) 

^0 



The terms up to one-loop (Dq, £>!, iVi) have been calculated [3[T0l[TTl[l2l[ia[Il],[T5l[ia[IIl[Il], 
[S], whereas only some parts of N2 are currently known [111 HO]- The inclusion of the NiDi/D^ 
term without N2 would hence be incomplete, and we have chosen to use only the lowest order cross 
section in the denominator and the 0{a^) term in the numerator, as done in [S|. 

Rewriting N and D to include the EW contributions yields the following expression for the 
leading terms. 



Apr = — = 



N a^No + a^Ni + a'^aNi + af 7V2 + 



Ni Ni 



D a'^DQ + alD^ + alDi+alaDi 



'Dn 



Dn 



a^ No 
allTo 



(6) 



where the incomplete 0{a1) part has been dropped. In the following we (re-)evaluate the three 
contributions on the r.h.s. of dfil). 





Figure 1: Electroweak and QCD Born diagrams 



Figure [T] contains all the tree level diagrams for the partonic subprocesses qq -^ tt and gg — >• tt 
(Higgs exchange is completely negligible). The squared terms \M.qq^g^tt\^ and \Mgg^u\'^ yield 



Dq of the LO cross section; the O(a^) terms arise from \A4qq-yj-^tt+ Mqq^z->ti\'^: which generate 
a purely-electroweak antisymmetric differential cross section, in the parton cms given by 
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In ApB © this leads to the term Nq- The complementary symmetric cross section provides the Dq 
term in the denominator, which does not contribute in the order under consideration. Interference 
of gg — > 7, Z ^ ti and qq ^ g ^ tt is zero because of the color structurqj. 

The 0{al) terms that contributes to N arise from four classes of partonic processes: qq — >■ ti, 
qq — )■ tig, qg — >■ tiq and qg — >■ tiq. In the first case the origin is the interference of QCD one-loop 
and Born amplitudes; the other processes correspond to real-particle emissions. All one-loop vertex 



q g , 



Figure 2: QCD box diagrams 

corrections and self-energies do not generate any asymmetric term, hence only the box diagrams 
(Figure[2]) are relevant. The box integrals are free of ultraviolet and collinear divergences, but they 
involve infrared singularities which are cancelled after adding the integrated real-gluon emission 
contribution qq -^ tig, shown in Figure [3l For the corresponding relevant gluon-radiation part 
only the interference of initial and final state gluon radiation has to be taken into account, yielding 
another antisymmetric cross section. The processes of real-quark radiation qg — > tiq and qg — > tiq 




Figure 3: Real emission of gluons at 0{a'l) 



^For qq — > tt there are also 0{a) W-mediated t-channcl diagrams with q = d,s, b, but they are strongly suppressed 
by the CKM matrix or by parton distributions [q = b). 



yield contributions to Aps which are numerically not important [5]. 

In order to analyze the electroweak 0{a'^a) terms, it is useful to separate the QED contributions 
involving photons from the weak contributions with Z bosons. In the QED sector we obtain the 
0{ala) contributions to N from three classes of partonic processes: qq -^ ti, qq — >■ tig and qq — >■ ti'y. 
The first case is the virtual-photon contribution, which can be obtained from the QCD analogue, 
namely the 0{a^^) interference of box and tree-level amplitudes, by substituting successively each 
one of the three internal gluons by a photon, as displayed in Figure H] 




q 9 t q 9 i. 
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Figure 4: Different ways of QED-QCD interference at 0{ala) 



The essential differences between the calculation of the 0{al) and of QED 0{a^a) terms are 
the coupling constants and the appearance of the SU{3) generators in the strong vertices. Summing 
over color in the final state and averaging in the initial state wc find for the virtual contributions 
to the antisymmetric cross section the following ratio. 



\M*t 



O (a'^a) ,asym 



2i?e(>l«(„)A^g(4))^^ 
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(8) 

that can be expressed in terms of two factors Fq^j^ and Fq^jj^ depending only on coupling constants 
and color traces, 



T7tt 
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(9b) 



Fqcd contains two different color structures and the result depends on (P = d'^^'-" dABC = -y, 
which arises from Tr{tH^t^) = jiif^^^ + d"^^^). Fq^^ instead depends on the charges of the 
incoming quarks (Qq) and of the top quark (Qt), together with n^j = 3 corresponding to Figure|H 



In a similar way, also the real-radiation processes qq 
can be evaluated starting from the result obtained for qq - 
successively each gluon by a photon, yielding the ratios 



> ttg and qq — > ^7 (Figures [S] and [5]) 
ttg in the QCD case and substituting 
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Fq^jj, Fq^jj and Fq^jj arc related to Fq^j^, Fq^j-, in the following way. 



(11) 
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This guarantees the cancellation of the IR singularities stemming from the virtual contributions. 

The 0{aj.a) antisymmetric term from qq — > tig comes from the interference of gg -H- g — >■ tig 
(Figure [3]) and qq ^^ ^ ^^ tig (Figure [5]). It can be obtained from the corresponding QCD result 
with the replacement of one gluon by a photon and the right couplings, as done in the case of 
qq -^ ti. The only difference is the number of gluons to be replaced: in the qq — >■ tig case they are 
only two instead of three as for the virtual photon contributions. 




Figure 5: Real gluon emission from photon exchange diagrams 

The ©(a^a) antisymmetric term from qq — >■ ti^ comes from the qq ^ g ^>- tij diagrams in 
Figure [6l and again it can be obtained from the corresponding QCD result for the gluon-radiation 
process qq -^ tig. Here we have a one-to-one relation between the QED and QCD diagrams. 

Finally, we can relate the QED contribution to the antisymmetric term A^i in ^ to the O(a^) 
QCD term A^i for a given quark species qq ^ ti+ X in the following way, 
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(13) 



Figure 6: Real photon emission from giuon exchange diagrams 



Now we consider the weak contribution to iVi. It can be depicted by the same diagrams as 
for qq — >■ ti and qq -^ tig in the QED case, but with the photon now substituted by a Z boson, 
involving massive box diagrams. The result cannot be expressed immediately in a simple factorized 
way. We performed the explicit calculation including also the contribution from real gluon radiation 
with numerical integration over the hard gluon part. 

Basically also Z-boson radiation, qq -^ ttZ , can contribute at the same order. As our calculation 
has shown, it yields only a tiny effect of 10~^ in ApB and thus may be safely neglected. The same 
applies to ud — > ttW^ as well as to Higgs-boson radiation. 

Weak one-loop contributions to the qqg and tig vertices induce also axialvector form factors, 
which however yield vanishing interference terms with the Born amplitude for the antisymmetric 
cross section at 0{ala) and are thus irrelevant. 



3 Numerical results 



The numerical analysis is based on the analytical evaluation of the required symmetric and an- 
tisymmetric parts of the parton cross sections and semi-numerical phase-space integration for the 
radiation processes with phase-space slicing, with support of FeynArts |21] and FormCalc [55j. This 
is done also for the QED subclass starting from the ti, tig and ti-j diagrams, for comparison with 
the QED result obtained from (fT3|) . showing perfect compatibility. 

We choose MRST2004QED parton distributions [23] for NLO calculations and MRST2001LO for 
LO [24], using thereby as(/z) of MRST2004QED also for the evaluation of the cross sections at LO 
(a similar strategy was employed in [5]). The same value fj, is used also for the factorization scale. 
The numerical results are presented with three different choices for the scale: /x = mt/2,mt,2mt. 
Other input parameters are taken from |25j . 

The results for the cross sections from the individual partonic channels and their sum, yielding 
the denominator of ApB , are listed in Table [T] The various antisymmetric terms entering the 
numerator of either of the two variants A"^ and A^p^ are collected in TableO and the corresponding 
contributions to the asymmetry in Table |4| 

As already mentioned, the QED part was obtained in two different ways based on the diagram- 
matic calculation and on the use of ([T3| : the weak part results exclusively from the diagrammatic 
calculation. The ratio of the total 0(0^0) + C(a^) and 0{al) contributions to the numerator TV 
of the asymmetry ^ gives an illustration of the impact of the electroweak relative to the QCD 
asymmetry. The values obtained numerically for /i = {mt/'2,mt,2mt) for the two definitions of 



a{pb) 


11 = ■mt/2 


tJ. = mt 


tJ. = 2mt 


uu 


6.245 


4.454 


3.355 


dd 


1.112 


0.777 


0.575 


ss 


1.37 X 10^2 


9.60 X 10-3 


0.706 X 10-2 


cc 


2.24 X 10-3 


1.69 X 10-3 


1.32 X 10-3 


99 


0.617 


0.378 


0.248 


PP 


7.990 


5.621 


4.187 



Table 1: Integrated cross sections at 0{al) from the various partonic channels 



ApB are 



TV*. 



R 



EW 



ATti 



= (0.190,0.220,0.254), 



RPP^^ = ^(^l^±£(f!I = (0.186,0.218,0.243), 



•■EW 



atPP 
C)(q3) 



(14) 



which are larger than the estimate of 0.09 given in [5]. This shows that the electroweak contribution 
provides a non-negligible fraction of the QCD-based antisymmetric cross section with the same 
overall sign, thus enlarging the Standard Model prediction for the asymmetry (the electroweak 



Oial 



contribution of ■ 



tt to the asymmetry is even bigger than the 0{a'l) contribution of 



dd^tt). 

The final result for the two definitions of ApB can be summarized as follows. 



Afg = (9.7, 8.9, 8.3)%, 



APP 



(6.4,5.9,5.4)%. 



(15) 



Figure [7] displays the theoretical prediction versus the experimental data. The prediction is 
almost inside the experimental la range for A"^ and inside the 2(T range for J^pg- It is important 
to note that the band indicating the scale variation of the prediction does not account for all the 
theoretical uncertainties. For example, the 0{a^) term in N is missing, and we did not include the 
0{a'l) part in D. Including the NLO term for the cross section in D would decrease the asymmetry 
by about 30%, which indicates the size of the NLO terms. In a conservative spirit one would 
consider this as an uncertainty from the incomplete NLO calculation (see also the discussion in [5]). 

We have performed our analysis also for applying two different types of cuts, one to the tt 
invariant mass and the other one to the rapidity: Mfi > 450 GeV and |Ay| > 1. With those cuts, 
experimental data have also been presented in [5]. The cross section values for these cuts at LO 
are given in Table [2] The various terms of the antisymmetric cross section contributing to N , as 



a(ph) 


A* 


= mt/2 


tJ. = mt 


11 = 2mt 


pp{Mti > 450 GoV) 




3.113 


2.148 


1.573 


pp(|A'y| > 1) 




1.846 


1.276 


0.937 



Table 2: Cross sections with cuts at 0{al) 
discussed above in the case without cuts, are now calculated for A"^ for both cases Mft > 450 GeV 



7 



and |Ay| > 1. The corresponding contributions to the asymmetry A"^ are the entries of Table [SJ 
The asymmetry with cuts is the total result, 

A%B{Mtt > 450 GeV) = (13.9, 12.8, 11.9)%, A^J^^dAj/l > 1) = (20.7, 19.1, 17.5)%. (16) 



A comparison of Table [5] with Table 3(a) shows that the ratio of the QCD contribution to the 



uu -^ tt and dd -^ tt subprocesses is larger with the M^f > 450 GeV cut, which leads to a slight 
increase of R*ew'- 

RswiMti > 450 GeV) = (0.200, 0.232, 0.266) R%w{\Ay\ > 1) = (0.191, 0.216, 0.246). (17) 

It is, however, not enough to improve the situation. 

Figure [8] displays the theoretical prediction versus data for Apg with cuts. The Standard 
Model prediction is inside the 2a range for the \Ay\ > 1 cut, but it is at the 3cr boundary for the 
invariant-mass cut M^j > 450 GeV. 

4 Conclusions 

Our detailed analysis of the electroweak contributions to the forward-backward asymmetry in tt 
production shows that they provide a non-negligible fraction of the QCD-induced asymmetry with 
the same overall sign, thus enlarging the Standard Model prediction for the asymmetry at the 
Tevatron. For high invariant masses, a 3tT deviation from the measured value still persists. The 
observed dependence of ApB on the invariant mass of ti could be an indication for the presence of 
new physics below the TeV scale; it is, however, difficult to interpret these deviations as long as the 
NLO QCD calculation for the asymmetry is not available. 
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Table 3: The various contributions to the antisymmetric cross section TV of Ap^ and A^j 
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Table 4: Individual and total contributions to A^pg and A^pg 
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Figure 7: Tlieory(blue) and experimental data (black=central value, orange^lu, yellow=2(T) 
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Table 5: Individual and total contributions to A^^^^Mti > 450 GeV) and A*^g{\Ay\ > 1) 
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Figure 8: Tlieory(blue) and experimental data (black=central value, orange^lc, yellow=2(7) 
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